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These notes provide an introduction to a number of those topics in Classical Mechanics that are 
^ . useful for field theory. 

1.1 Lagrange's Equations and The Action Principle 

In analysing a physical system, the principal goal is to study its time evolution. If the state of a 
^ ■ system is known at time t=0, we wish to know how it will appear at subsequent times. For Newton, 
this meant that the positions rj(t) and velocities Vj(t) = ri(t) of a collection of particles i = 1, N 
would evolve in time according to his famous equation F = ma. Since a(t), the acceleration, is the 
time derivative of the velocity, this is a second order equation and hence one needs to have both 
the initial positions and velocities of all particles to see what happens to a system as times elapses. 

For the particular case of a conservative force field, the force F can be written as the gradient 
of a time independent and velocity independent scalar function V, the potential, so that Newton's 
equation reduces to 

mi = -W. (1) 

Upon multiplying this equation by r , this equation can be integrated once with respect to time to 
yield the result 

E = imr 2 + V = T + V (2) 

where E, the energy, is a conserved quantity. (Remarkably, Newton himself didn't use conserved 
quantities in his analysis of dynamical systems; he worked entirely from the equations of motion.) 
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In passing to the Lagrangian formalism, one makes use of generalized coordinates qi(t). The 
space of these coordinates is generally called "configuration space". It is easily shown that if 
r, = Ti(qi...q n ), then Newton's equation can be rewritten in the form 

d dL dL n 
dt dqi dqi 

where L = L(q, q), the Lagrangian, is T — V. (Including explicit time dependence in L is trivial.) 
In order to show that eqs. (1) and (3) are equivalent, one first notes that eq. (1) can be written 

as 

d dT dV () 

dt dfj drj 



Then, since = implies that ^ = we see that 



dT _ dTdf^ _ dTd^ 
dqi dfj dqi drj dq t 



(5) 



Consequently, we have 



which by eq. (4) becomes 



d dT _ f d dT\ drj dT dfj 
dt dqi \ dt dfj ) dqi dfj dqi 



dV dr dT d 

3 + %- = Jf-(T-V). (7) 



dr j dqi dqi dq, 

This is equivalent to eq. (3). 

The solution of eq. (3) can be seen to extremize the action integral 

S= I"* L{ qi {t)Ai{t))dt. (8) 

To show this, let 

qi(t) = qcu(t) + e8 qi (t). (9) 

with 5qi(t\) = 5qi(t 2 ) = 0. (The subscript CL refers to the "classical" trajectory.) We then have 
the action S depending on the parameter e. The requirement that ^ vanishes when e equals zero 
leads to 

r* WfaW faW) _ d9L(M,Jt)\ {t)dt = (10) 

Jti V %W dt dQi(t) J 

If this were to vanish for arbitrary 5q(t), then q Ch (t) would have to satisfy the Lagrange eq. (3). 

When dealing with this second order differential equation, one normally specifies the initial 
values of qi and in the variational approach, the initial and final values of are fixed. In both 
cases, two boundary conditions are applied for each value of i. 



In principle, the Lagrangian could depend on derivatives of qi(t) beyond the first. If, say, L = 
L(qi(t), qi(t), (ji(t)), then the integral Ldt would have an extremum on the trajectory satisfying 
the fourth order equation 

OL _ d OL d 2 dL 

dqi dt d<ji dt 2 0q\ 

provided that the initial and final values of qi and are specified. This extension of Lagrange's 
equations was first considered by Ostrogradsky. 



1.2 Hamilton's Equations and Poisson Brackets 

In the Lagrangian approach to dynamics, the only independent variables are the N components of g«. 
If L depends only on q^ and <ji, then the time evolution of is dictated by N second order coupled 
ordinary differential equations. It is possible in this case to introduce a further iV independent 
variables pi (the canonical momenta) and through a Legendre transformation arrive at a set of 2N 
first order ordinary differential equations in the 2N dimensional "phase space" of the g's and p's 
that determines their time evolution. 
One begins by defining 

Pi ~ 777- (12) 

and then introduces a Hamiltonian H using the Legendre transformation 

H(qi,Pi) = ^Pitji- L(q i ,q i ). (13) 

i 

(We henceforth will use the Einstein summation convention for repeated indices.) It is understood 
for the present that we are able to eliminate dependence on the by expressing them in terms of 
the Pi. It follows that, H does not have dependence on q^ as W- vanishes on account of the way in 
which pi and H have been defined in eqs. (12) and (13). Variation of H can be expressed in two 
ways; assuming that H has no explicit time dependence, first we see that 

OH OH 

dH = -r^dqi + -z-dpi (14) 
Oqi Opi 

and then from the definition of H 

OL OL 

dH = tjidpi + pidqi - — dq t - T^dtji. (15) 

Oqi 0q { 

In eq. (11) we can replace ^ by p (from the definition of p) and ^ by p (from Lagrange's equations). 
Matching coefficients of dp and dq in eqs. (14) and (15) then yield the 2N Hamilton's equations 



dH . . 

Pi = 17 

We first note that if L, and hence H, has explicit time dependence, then 

dH _dH . 9//. <9# 

rft 9g« ^ dpi ^ % dt 

On account of the Hamilton's equations of motion, H is a constant in time provided it has no 
explicit time dependence. If 

L = ^Mijqtqj - V( qi ), (19) 

then it follows that H = T + V, numerically equal to the total energy E when the equations of 
motion are satisfied, with T = \M i) jq i qj being the kinetic energy. 

Hamilton's equation can in fact be derived by requiring that the action integral 

*t 2 



[\p i q i -H(q i ,p i ))dt (20) 
Ju 



be an extremum under independent variations of both ^ and Pi, provided that qi is fixed at the 
initial and final times; pi need not be fixed in this way. 

It is convenient to define the Poisson bracket of two dynamical variables A(q,p) and B(q,p) to 

be 

, , n1 fdAdB dBdA\ 

i A ^}pB= o-t--t-t- ■ (21) 



\dqidpi dqidpi 

(The subscript "P_B" will henceforth be written explicitly only if there is an ambiguity.) 
It follows from Hamilton's equations that 

dA(q(t),p(t)) _ dA OA 
dt dqi dpi 

can be written 

£ = (23) 

This equation allows us to reduce much of dynamics to an algebraic problem as the Poisson brackets 
satisfy some general relations that follow from their definitions. We note that: 

* - Uiilj} = = {Pi,Pj} ( 24 ) 

ii-{q i ,Pj} = 5 iij (25) 

Hi - {A,B} = -{B,A} (26) 

iv-{A + B, C} = {A, C} + {B, C} (27) 

v — {A, BC} = B {A, C} + {A, B} C (28) 

vi - {A, {B, C}} + {B, {C, A}} + {C, {A, B}} = 0. (29) 



Only property vi, the Jacobi identity, is non-trivial to prove. From this identity, it follows that if A 
and B both have vanishing Poisson brackets with the Hamiltonian, then so does {A, B}. 

The canonical formalism can also be developed to deal with the situation in which the Lagrangian 
depends on q as well as q and q. In this case, we make the definitions v = q, p = ^ — ^ (jfj^j an d 
7r = eq. (11) can then be expressed as a set of first order equations q = v = p = — ^ 
and 7r = — ^ where the Hamiltonian is a function of the independent variables q, v, p and ir given 
by H(q,v;p,ir) = pq + ixv - L(q,q,q). 



1.3 Canonical Transformations and the Hamilton- Jacobi Equation 

A change of variables from (qi,Pi) to (Q i: Pi) is said to be "canonical" if it leaves the form of 
Hamilton's equations unchanged. Thus if (qi,Pi) satisfy Hamilton's eqs. (eqs. (16) and (17)), there 
exists a function K(Qi, Pi) such that 

dK 

Q, = ^ (30) 
^ 3K (3D 



The action S = f**(PiQi — K) is extremized by the solution to these equations, just as the 
action of eq. (20) leads to Hamilton's equations for qi and pi. As the variations of (q,p) and (Q, P) 
vanish at the end points ti and t 2 , these two integrands can differ at most by the time derivative of 
some function F, so that 

p t q t -H = PiQ t -K + ^. (32) 
The relation F can be function of only 2N of the AN variables (qi,Pi, Qi, Pi); we initially examine 

F = F 1 (q l ,Q i ,t). (33) 

(We now consider the possibility that H has explicit time dependence.) As 

dF dF 1 . 3Fii dFi 

-dT = ^q- qi+ dQ^ + ~dt (34) 



together eqs. (32) and (34) imply that 



If now we were to set 



K = H +1 ± (37) 



F = Fi(q„P i )-Q,P i (38) 



then it follows in a similar fashion that 

dF 2 



% 

Q* = U (40) 



K = H + ^. (41) 



Next, if we take 



we obtain 



di 
~dt 

F = F 3 (p t ,Q l ) + q tPl (42) 



ft = (43) 



- (44) 
K = H + ^.. (45) 



Lastly, there is the case in which 



this results in 



dp, 
dF 
dQ 

m 

~dt 

F = F 4 + q iPi - Q^; (46) 

* = -<sr (47) 
* = i (48) 

K = H + § . (49) 
The transformation generated by F 2 is the identity if 

F 2 = g^; (50) 

an infinitesmal canonical transformation takes place if 

F 2 = q i P i + eG(q i: P i: t) (51) 

where e is a small parameter. In this case eqs. (39) and (40) result in 

dG 

5pi = Pi-pi = -e— (52) 
% 

dG 

5qi = Qi- qi = e—. (53) 

Upon taking e = 8t, a small increment in time, and identifying G with the Hamiltonian H, these 
equations are seen to be identical to Hamilton's equations, so that the time evolution of a system 



can be viewed as sequence of canonical transformations. Similarly, a translation in a direction i by 
an amount e is generated by the momentum pi (that is, we just take G = Pi ~ Pi), and a rotation 
about an axis i through an angle e is generated by the angular momentum 



Li = ei jk qjp k . (54) 

The fundamental Poisson brackets lead to the interesting result that 

{Li, Lj} = eijkLk- (55) 

The notion of a canonical transformation can also be used to reduce a dynamical problem to that 
of solving a single partial differential equation, the Hamilton- Jacobi equation. The Hamilton- Jacobi 
function is a transformation function F 2 (qi, Pi,t), often denoted by S, that results in a vanishing 
new Hamiltonian K. Eq. (41) then reduces to 

In this equation, we have used eq. (39) to eliminate Pi in the argument of H(q,i,pi) in favour of 

The momenta Pi are necessarily time independent as the new Hamitonian K vanishes; these 
momenta are identified as being the constants of integration that arise in the course of solving for 
S in eq. (56). Calling these constants cti, we see that eq. (40) implies an additional set of constants 
given by 

* = S- < 57 » 

Once we have solved eq. (56) for S we can then use eq. (57) to obtain 

qi = qi(ati,Pi,t). (58) 

Thus the constants a and /3 can be identified with the boundary conditions on the initial configu- 
ration of the system that are needed in the course of solving for its dynamical evolution. 
A formal result is that 

dS( qi (t),ai,t) _ OS . OS 

It " dq~ qi dt (59) 

which by eqs. (40) and (56) yields 

S= J(piqi-H)dt. (60) 

This equation shows that the action S appearing in eq. (20) is in fact identical to the Hamilton- 
Jacobi function. This doesn't make it possible to find the Hamilton- Jacobi function though, as the 
integral of eq. (60) can only be evaluated once g«(t) and pi{t) are known. 



As a simple example, let us consider a one dimensional harmonic oscillator in which both the 
mass and angular frequency are scaled to one. In this case, the kinetic and potential energies are 
given by 



so that with L = T — V and H given by eq. (13), we have 

H=\tf + q% 
If the Hamilton- Jacobi function S is taken to have the form 



S(q, a, t) = W(q, a) — at 



(61) 
(62) 

(63) 
(64) 



then eq. (56) in this case reduces to 



dW 
dq 



q 



a. 



(65) 



(The constant a is evidently the energy E of the oscillator.) As a result of this equation, we find 
that 

(66) 



S — J dq\/2a — q 2 — at. 



Now using eq. (57), we arrive at 



dq 



V^a - q 2 



-t 



which can be integrated to yield 

q = V2asin(t + (3). 
Hamilton himself recognized that in the case where 

H=^ + V(q) 



(67) 
(68) 

(69) 



the Hamilton- Jacobi equation is related to the wave equation. If a wave travels in empty space with 
velocity c and is in a medium with refractive index fj,(q) , then the wave equation is 



H d ip 



o. 



(70) 



Restricting our attention to a single frequency u> = kc, we examine solutions to eq. (70) that are of 
the form 

^j(q,t) = M<iy {kf{q) ' UJt) - (71) 



For large values of k (in other words, short wave lengths), eq. (70) then reduces to 



(V/) 2 = fi 2 . (72) 

If one were to take the form of the Hamilton- Jacobi function S associated with the Hamiltonian of 
eq. (69) to be 

S = akf(q) - Et (73) 
then the Hamilton- Jacobi eq. (56) with this Hamiltonian becomes 

J-(aWf(q)) 2 + V(q) = E (74) 
which is the same as eq. (72) provided 

2m(E - V) 

" = aW ■ (75) 
This identification means that the original wave eq. (70) would become 



The constant a can now be identified with h so that eq. (76) becomes the time independent 
Schrodinger equation. 



1.4 The Dirac Constraint Formalism and Gauge Invariance 

In some cases it turns out that it is convenient to consider Lagrangians in which not all of the 
variables in configuration space are independent. (This may be to ensure that symmetries present 
in the system being considered are readily apparent.) In this case the equations of motion do not 
uniquely determine the evolution of the system as time elapses. The Dirac constraint formalism 
provides a systematic way of treating such systems. 

To appreciate the sort of difficulties that can arise in extracting information about the time 
evolution of such systems, consider a Lagrangian L(qi(t) , (i = 1...N) with ?j(0) and <&(0) 

being prescribed initial conditions. If we now examine a Taylor series expansion of qiit) about 
t — 0, then 

qi(t) = fc(0) + tqm + h 2 q t (0) + (77) 

with the boundary conditions providing the first two terms of this expansion. In the third term, the 
value of (ji(0) is in principle fixed in terms of g;(0) and ^(0) by Lagrange's eq. (3). This is because 
eq. (3) can be expessed in the form 



allowing us to solve for q\(0) provided the so-called Hessian matrix 

d 2 L 

can be inverted. (Taking the time derivative of the Lagrange equation allows us in principle to 
find higher derivatives of % at time t — 0.) Systems in which this inversion is not possible contain 
constraints and do not have the values of (ji(0) (and consequently of fixed in terms of the 

boundary conditions. 

Inability to invert the Hessian matrix amounts to being unable to solve eq. (12) so that the 
generalized velocities are expressed in terms of the canonical momenta Pi without the imposition 
of a set of "primary" constraints 

K t (q,p)=0. (80) 

The number of these primary constraints is N — R where R is the rank of the Hessian matrix. 

Once these constrants are satisfied, it becomes possible to define the canonical Hamiltonian He 
as in eq. (13); however, in order that eq. (80) is satisfied we consider the "total" Hamiltonian 

H T = H c + u i 7r i (81) 

where the fields set of Lagrange multiplier fields. The equations of motion for these fields 

ensure that the primary constraints are satified. 

However, for consistency it is also necessary to have 

^ « (82) 
dt 

where ~ denotes a "weak" equality, one which holds when the constraints themselves are satisfied. 
Using eq. (23), this means that 

{n t , H T } « 0. (83) 

This condition may be automaticly satisfied. It may also be solved for some or all of the Lagrange 
multiplier Ui in eq. (81). Lastly, it may require imposition of additional constraints 

<Ti(q,p)=0, (84) 

where these Oi are known as "secondary" constraints. The process applied to the primary constraints 
is now repeated; we consider 

KHt}«0 (85) 

and check to see if the time derivative of all of constraints, both primary and secondary, vanishes 
weakly. Further tercery constraints may arise (though this does not happen often). One continues 
adding constraints until all constraints have a vanishing Poisson bracket with the Hamiltonian when 
the constraints themselves are satisfied. 



Having obtained all constraints (ipi), it is convenient to divide them into first class constraints 
((pi) and second class constraints A first class constraint has a weakly vanishing Poisson bracket 
with all other constraints so that 

{0,,^}~O; (86) 

a second class constraint is one that is not first class. (Separating the first class constraints from 
the second class constraints can be difficult.) This distinction between the two classes of constraints 
is important, as the first class constraints are associated with arbitrariness in the time development 
of the system while second class constraints can be used to eliminate degrees of freedom in phase 
space. To see this, we consider the extended Hamiltonian H E in the form 

H E = H C + Ci<f>i + di9 t (87) 

with the Lagrange multipliers c and d associated with the first and second class constraints respec- 
tively. The requirement 

{0,, H E } » (88) 

is satisfed for all values of q and di as first class constraints have a weakly vanishing Poisson bracket 
with both first and second class constraints. However, when we consider 

{6i, H E } » (89) 

we find a set of conditions on the Lagrange multipliers df, otherwise d; L would be first class. If now 
we apply eq. (23) to examine how a dynamical variable A evolves in time, we find that since the 
Lagrange multipliers q appearing in He are arbitrary, ^ is not fixed. 

Elimination of superfluous degees of freedom using second class constraints involves use of Dirac 
brackets in place of Poisson brackets. This entails first defining the antisymmetric matrix 

d,.,^ {0,.0 s }. (90) 

This matrix is invertible; otherwise not all of the constraints di would be second class. The number of 
second class constraints is consequently even because the dimension of any invertible antisymmetric 
matrix must be even. We now define the Dirac bracket of two dynamical variables A and B to be 

{A, B}* = {A, B} - {A, OA d~j ft, B} . (91) 

These brackets have a number of useful properties. First of all, we note that for any second class 
constraint 0, we have 

{6i,By = 0. (92) 

Next, if F is a first class dynamical variable (that is, one whose Poisson bracket with any constraint 
ipi weakly vanishes so that {F, tpi} 0), then it follows immediately that 

{F,B}*{F,B}*. (93) 



The Hamiltonian is such a first class quantity, as can be seen from eq. (23) and the consistency 
condition that all constraints have a vanishing time derivative. As a result of eq. (93) then, we can 
replace eq. (23) by 

HA 

^~{A,H T y ; (94) 

this equation need hold only weakly, that is, on the constraint surface. But now eq. (92) tells us 
that in Ht we can set the second class constraints 6i equal to zero. It is consequently not necessary 
to determine the value of the Lagrange multipliers di appearing in eq. (87). Thus by use of Dirac 
brackets, we have been able to eliminate some of the degrees of freedom in phase space through 
setting the second class constraints equal to zero. 

The Dirac brackets can be shown to satisfy eqs. (26-29); proving the Jacobi identity (29) for 
the Dirac bracket is a tedious exercise. 

It is still necessary to deal with the arbitrariness noted above that would follow from the possible 
presence of first class constraints. Rather than simply choosing a value for the Lagrange multipliers 
Ci appearing in eq. (87), we introduce extra conditions, known as gauge conditions, in phase space 

7i(fc,Pi) = (95) 

with one gauge condition for each first class constraint. (The reason for calling these extra conditions 
"gauge conditions" becomes apparent later when the Dirac constraint approach is used to analyze 
the canonical structure of the electromagnetic field.) These gauge conditions are arbitrary functions 
of qi and pf, they are independent of t and the velocities qi and pi. We also require that these gauge 
conditions, together with the first class constraints (pi, form a set of second class constraints. In 
other words, we impose the restriction that if 

Di,j~{0i,7j} (96) 

then 

deW j£ 0. (97) 

As a result, the first class constraints, the second class constraints and the gauge conditions together 
form a set of constraints which can be treated as the second class constraints were by themselves. If 
9j denotes the set {0j,^,7j}, then Dirac brackets can be defined as in eq. (91), with 6j replacing 
9i. 

We now examine a simple system that serves to illustrate how Dirac brackets can be used. 
Consider a system in which a particle is constrained to move in a circle of unit radius with constant 
speed. Obviously, the easiest approach would be to use the angle 8(t) which gives the position 
of the particle as the only variable in phase space. Instead though, suppose we perversely use 
the Cartesian coordinates x(t) and y(t) to describe the position of the particle and in this two 
dimensional configuration space consider the Lagrangian 

L = i(i 2 + y 2 - x 2 - y 2 ) + \ 1 (x-y) + A 2 (y + x) . (98) 



The fields Ai and A2 are Lagrange multiplier fields whose equations of motion ensure that the 
particle moves in a circular path at constant speed; their presence means that configuration space 
is now four dimensional. 

The momenta associated with x, y, \± and A2 are 

p x = x + Ai (99) 

p y = y + A 2 (100) 

Px 1 =0 = p X2 . (101) 

The last equation constitutes a pair of primary constraints as one cannot solve for either Ai or A2 
in terms of p\ ± and p\ 2 . 

However, if we use these primary constraints, we can define the canonical Hamiltonian 

H c = Pxi + p y y + PA1A1 + p\ 2 \ 2 - L (102) 

and find that with the Lagrangian of eq. (94) 

H c = \ [{p x - Ai) 2 + {p y - A 2 ) 2 + x 2 + y 2 ] + X lV - \ 2 x. (103) 

Applying eq. (83) with the primary constraints of eq. (101) leads to a pair of secondary 
constraints 

Px - Ai - y = = p y - A 2 + x. (104) 

Together, eqs. (83) and (104) make up a set of four second class constraints, as {Aj,^ } = 
It is quite easy to compute the matrix d it j of eq. (90) associated with these four constraints; once 
that is done, it is possible to define the Dirac bracket of eq. (91) and we finally end up with the 
fundamental Dirac brackets 

{x,p x }* = 1 = {y,p y }* (105) 

with all other fundamental Dirac brackets vanishing. The Hamiltonian of eq. (103) itself, upon 
treating the constraints as strong equations, reduces to simply 

H c = Pxy - P y x. (106) 

With this Hamiltonian, we find that 

x ~ {x, H}* = y (107) 

and 

y^{y,H}* = -x (108) 



as anticipated. 



The simple system just considered does not contain any first class constraints. Before dealing 
with a particular situation in which first class constraints do occur, let us examine the nature of 
the gauge transformations arising due to the presence of first class constraints. Suppose we restrict 
ourselves to the situation in which there are only first class constraints, and these are either primary 
(nli) or secondary (vr2j), with 7r2j arising from the requirement that vanish when the primary 
constraints are zero. Inclusion of tercery constraints and beyond, as well as second class constraints, 
is quite easy. 

The arbitrariness associated with the presence of first class constraints means that trajectories 
in phase space are not unique; trajectories given by qi (t) and Pi (t) as well as by q^t) + and 
Pi(t) + /3i(t) can both be physical. We now take the generator that relates these two trajectories to 
be a function G whose form is given by 

G = e(t)G + e(i)Gi (109) 

where e is an infinitesmal parameter that may contain time dependence. (If there were any tercery 
constraints present, then a term dependent on e(t) would also contribute to G.) As this function G 
is a generator, we have the equations 

f)G f)G 

^) = {^G} = e(t)^ + e(t)^ (110) 

and 

r)G r)G 

m = { Pl ,G} = - e (t)-^-e(t)^. (Ill) 

From eq. (110) we see that 

a^e— + e — ,ff T +e— + e — ,F r . 112 
dpi \ dpi J dpi y dpi J 

(We take the Hamitonian to only incorporate the primary first class constraints in this discussion; 
weak inequalities hold on the constraint surface 7rlj = and not necessarly when 7r2j = 0.) The 
weak version of eq. (16) also holds for both qi and qi + af, we hence find that 

. . dH T (q i + a i ,p i + fa) 

% + ^ • ( 113 ) 

dpi 

To lowest order in ai and eq. (113) leads to 

8 2 H T a d 2 H T . . 

opidqj opiOpj 

We can equate dj found from eqs. (110) and (114) and then use eqs. (110) and (111) to eliminte 
ctj and fa to obtain 

^- [e {Go, H T } + e(G + {G ± , H T }) + ed] » 0. (115) 



Working with the equation for that is analogous to eq. (114), it is also possible to find that 

^- [e {Go, H T } + e(G + {G u H T }) + eG,} « 0. (116) 

Together, eqs. (115) and(116) can be solved by taking 

Gi = vrlj, (117) 

and then if 7r2j is found by considering {nli, H T }, 

Go = (-tt2 < + AtfTrlj) (118) 

where Ay is finally determined by the requirement that 

{G ,#t}~0. (119) 

We see that there is a separate generator G and associated gauge function e for each of the first class 
primary constraints 7rlj. This way of finding the generators of gauge transformations will prove to 
be quite useful when we come to analyze field theories containing first class constraints. 

An alternate approach to finding the gauge invariance of the classical action involves considering 
the "extended action" ^ 

Se= r dtipiqi- H E ) (120) 



Jti 

with 

H E = H C + U ai (j> ai + U ai e ai (121) 

where (f> ai is an i th generation first class constraint and 6 ai is an i th generation second class constraint. 
We know that 

{0 a , <M = G> c + ClXh + (122) 

{4>M = CiJ> + C! e fl fi (123) 

and that 

{H C , <p a } = v a b <f> b + vfe a e p + vy a (124) 

{Hc,e a } = v^ b + v^e p . (125) 

(The form of G" 6 and V@ is constrained as {<p a , <pb} and {H c , 7 a } must be first class.) If now and 
Pi undergo the "gauge transformations" 

4. = «-(*)^ ta-^M^ (126) 

while 

5 £ C/ a = e a + t/ c e fc G 6 a c + C/ Q e 6 G 6 a a - e b V b a (127) 



5 e U a = U c e b (C£ + C^e^j - e b (y b a + V^O^j + U?e h C^ (128) 

then 

5 e ( Pl q t ~H E ) = j t (-e a 4> a + Pie a ^A - (129) 

Consequently, if e a (ti) = e a (t 2 ) = 0, the action of eq. (120) is unaltered by a gauge transformation. 
Eqs. (126-128) can be recast into the form 

5F = {F,fi a <j) a } (130) 

where [i a now depends not only on t but also p iy q iy U a and U a . In this case 



su* = f/y (<- + cgfy) - // (V + vr%) + c/y eg, , (132) 

D d +(u^ + U^ + ^+(^ + U«° + ..) (133) 



where 

Dt dt ' dU a ' " dt/ a " 7 ' V~ du ° ' " <9f/ Q 
is the total time derivative, exclusive of dependency on time through and pj. It then follows that 



5 (pitji - H E ) = ^ 



d 

dpi 



(134) 



The transformation that leaves the action of eq. (4) invariant can be deduced by considering 
the invar iance of t 

St= ^ dt (p t q t -H c - U a ^ ai - U ai 6 ai ) (135) 

where only primary constraints occur in the sum appearing in St- As Se of eq. (120) reduces to St 
of eq. (135) upon setting U a ' 2 . . . U am , U a ' 2 . . . U an equal to zero, the invariance of Se is determined 
by choosing the "gauge" in which U a2 = . . . = U a " = = 5U a2 . . . = 5U an and iteratively solving 
for fi ai . . . fi am using eqs. (131,132). The number of gauge functions fi a then equals the number of 
primary first class constraints and their time derivatives occur if there are secondary, tercery etc. 
first class constraints. 

The procedure used to fix fi ai involves equating the coefficients of the constraints arising in eq. 
(131) to zero after setting U a \ U a % 5U a \ 5U a > (i = 2 . . . N) to zero. If % = 1...N (i.e. there 
are N generations of constraints) then we first take fi aN to be solely dependent on t, and then 
determine /x ai . . . n aN - 1 are found in terms of fi aN . These coefficients are uniquely determined by 
this procedure provided there is the same number of constraints occurring in each generation. This 
does not always happen; in a model in which a scalar field is on a curved surface in 1 + 1 dimensions 
there are two generations of constraints with more primary than secondary constraints and so the 
coefficients /z ai are not fixed by fi a ' 2 . 



1.5 Grassmann Variables and their Canonical Formalism 



Grassmann variables have the distinctive property that they anticommute under multiplication. 
That is, if 9\ and 2 are Grassmann variables, then 

0102 = -0 2 0i (136) 

As a result, we cannot "count" with Grassmann variables, and they are not ordered; no meaning 
can be given to an inequality such as 0i > 2 . (For Hermitian conjugation we have (0102^ = 9\9[.) 

Any Taylor series expansion of a function of Grassmann variables in powers of these variables 
must terminate after a finite number of terms since if is Grassmann, then 2 = on account of 
eq. (136). Thus a function F(x; 9 1: 2 , ■■■0 n ) can be expanded 

F(x;9 u ...,0 n ) = f (x) + /i(x)j0j + f 2 (x) ili2 9 h 9 i2 + ... + /„(x) ili2 ... in il i2 ...0 in (137) 

where the functions fm{x)i 1 ___i rn are a ^ antisymmetric in the indices %\...% m . 

One can introduce a calculus of Grassmann variables. Differentiation is defined so that 

w?* = 5ij - (138) 

It is also necessary when taking derivatives to take into account the fact that Grassmann variables 
anticommute, so for example ^0102 = —9\ and F = — de f ddi F. We also have j- t F{9{t)) = 
9(t)F\9{t)). 

Integration is defined to be the same as differentiation; we take 

fM S Jnn.W=4gJ^F (139) 

so that, for example, J d9id9 2 (9i9 2 ) — — 1. It also follows that a shift of variables leaves a Grassmann 
integral unaltered so that J d99 = J d9(9 + 0) where <fi is also Grassmann. 
Under a change of variables 4>i — a ij@ji it is evident from eq. (119) that 

J d n 9F{9) = J d n (pdet(a ij )F((f)) (140) 

and hence d n 9 = d n (pdet(aij). If and <fi were ordinary variables rather than Grassmann, then the 
usual rules of calculus give d n 9 = d n 4>det~ 1 (a i j). 

If x is a real vector, then the standard Gaussian integral d\e~ aX2 = easily leads to 

n 

d n xexp (-x T Ax) = ^=L= (141) 

for a (symmetric) matrix A. We also find that if z is a complex vector and H is an Hermitian 
matrix, then 

d n zd n z*exp (-z T *Hz) = (142) 



f 

J — c 



It is also possible to deal with Gaussian integrals over Grassmann variables. If 9 is an n com- 
ponent Grassmann vector with n even, and B an n by n antisymmetric matrix, then eq. (139) can 
be used to show that 

J d n 9exp (^-h) T B6^j = VdrtB. (143) 

Proving this involves making a transformation on B so that it has a sequence of 2 by 2 blocks along 
its diagonal with zeros everywhere else. (See eq. (2.166) below.) Furthermore, we can show that if 
9 and <p are independent Grassmann vectors of dimension n, then 

J d n <pd n 9exp (-fHO) = detH (144) 

for a matrix H . 

There is also a canonical formalism associated with Grassmann quantities that are dynamical 
variables. A Lagrangian in this case is a function of both ordinary variables q^t) and Grassmann 
variables 9i(t). Variation of the action 

S= ^ dtL{ qi {t)Ai{t)-A{t)Mt)) (145) 
Jti 

leads to the expected equations; eq. (3) as well as analogous equation for 0j. If we now define a 
Grassmann momentum conjugate to 0j , 

* = § (146) 

as well as a Hamiltonian H 

H(qi,Pi', Oi, 7Tj) = qiPi + Oi'Ki — L (147) 

then variation of the action ( 

S = jT 2 dt [q iPi + km - H) (148) 

leads to the canonical equations of motion of eqs. (16) and (17) as well as the equations 
and 

dH , . 

*' = ~w {m 

The unexpected minus sign in eq. (149) comes about as 9^6^ = —5iTi9i. 

In order to retain eq. (23), we define Poisson brackets to be now given by 

SA B\ -(——-—— \ - (—— 

\dpidqi dqidpi) \dQi d-K i d9 i ) ' 

The order in which terms appear in eq. (151) is important, as A and B may be Grassmann, as well 
as 9 and n. The ordering of Grassmann quantities can also lead to extra minus signs appearing in 
eqs. (26-29) when considering these generalized Poisson brackets. 



1.6 Special Relativity 

Any physical event is taken to occur at a specific time and place, as described by Cartesian co- 
ordinates (t,x,y,z). These coordinates are those used by a specific observer; some other observer 
moving with relative velocity v would use different coordinates which we take to be (t', x', y' ', z'). If 
the two observers are coincident at t — t' — with their axes are aligned at that time, and v is in 
the direction of the x and x' axes, then classically we have the so-called Galilean transformation 

x = x — vt, (152) 

v = y, (153) 

z = z, (154) 

and 

t' = t (155) 

where v is the magnitude of v. This transformation is consistent with everyday experience; in 
particular, it implies that if u — ^ and u — then 

u = u — v . (156) 

However, as we shall show later on, the Maxwell equations which describe electricity and magnetism 
predict the existence of waves travelling with velocity c. If c is to be universal constant, then this 
velocity must be the same according to all observers and eq. (156) breaks down. Einstein pointed 
out that this implies the relationship between the primed and unprimed coordinates used to label an 
event can no longer be given by eqs. (152-155); in particular the spatial and temporal coordinates 
must mix. (The actual transformation was discovered earlier by Voigt, Larmor and Lorentz, but 
they only viewed this as a useful invariance of the Maxwell equations.) 

If the clocks and coordinate axes coincide in the primed and unprimed frames at t — t' — and 
the relative velocity v is again along the x and x axes, then we assume the linear transformation 



and 



First of all, the point x 



x = aux + a^t, 

y = y, 

i 

Z = 

t = a 2 \X + a 2 it- 
= coincides with x = vt and so eq. 

a 12 = -va xl . 



(157) 

(158) 
(159) 

(160) 

(157) implies that 

(161) 



Since c is to be a universal constant, a wave front emitted along the positive x axis at time t — 
will be at x = ct and x = ct' in the unprimed and primed frames respectively. Equations (157) 
and (160) together with eq. (161) then imply that 

ct = a u (c-v)t (162) 

and 

t = (a 2 ic + a 2 2)t. (163) 
So also, a wave front along the negative x axis leads to 

-ct =a u (-c-v)t (164) 

and 

t = (-a 2 ic + a 22 )t. (165) 

From eqs. (162-165) it follows that 

x —an(x — vt), (166) 

t =a u (~x + t). (167) 
cr 

But because of the symmetry between our two reference frames, we must also have 

x = an(x + vt ) (168) 



and 



together eqs. (166-169) lead to 



t = a u (^x +t'); (169) 



c 2 



a n = * = 7 . (170) 



These transformations are generally called the "boost transformations" . 

f ^d u' = g 



It immediately follows that if u = ^ and u' = ^ 



u = —^L (172) 

u, = — -. — r . (173) 

7(1-^) 

These transformtion equations are consistent with the speed c being a universal constant. We shall 
henceforth scale c to one. 



If two events occur at the same position x and at times t\ and t 2 in one frame, then according 
to eq. (166), 

t 2 - t[ = 7 (t 2 - h) >t 2 - h; (174) 

this is so-called "time dilation". Similarly, if two events occur at the same time t' but at two 
different positions x\ and x' 2 (with the same y and z coordinates), then by eq. (168) we have 

i i x 2 — X\ 
x 2 — x 1 = < x 2 — X\ (175) 

7 

which is known as "length contraction" . 

These two phenomena lead to the so-called "twin paradox" . Suppose a rocket travels from earth 
to a star a constant distance D from earth then immediately returns, always moving with speed 
v. According to an observer on earth, the time taken for this trip is simply However, because 
an observer on the rocket is moving with respect to the earth and the star, he sees the distance 
between the earth and star to be just D\/l — v 2 on account of eq. (175), and so that according 
to this observer the time for the round trip is 2D ^^- V ' 2 . But by eq. (166), since the departure and 
return of the rocket both take place at the same place when viewed from the earth-bound frame, 
this time in the rocket's frame corresponds to a time ^ 2Dv/ J~" 2 ^ \/l — v 2 = (1 — v 2 ) ^ ^ ^ in the 
earth's frame. This inequality at first appears to be paradoxical; it can be accounted for though by 
noting that there is a sequence of events on earth during a time 2vD that are never simultaneous 
with either the out-bound observer (they are in his future) or the in-bound observer (they are in 
his past). 

To see this in more detail, let us call the earth-bound, out-bound and in-bound observers E, O 
and / respectively, so that according to eq. (167) 

to = l{-vx E + t E ) (176) 

and 

ti = i(vx E + t E ) (177) 

if the three frames used by these observers coincide at the instant the rocket initially leaves earth. 
Let us now consider three events; event L is the landing of the rocket on the star and its subsequent 
take-off towards earth, event A is an event on earth that is simultaneous in the frame O to L, and 
event B is an event on earth that is simultaneous in the frame I to L. Event L has coordinates 
(t E ,x E ) = (D,^). Furthermore, by eq. (176), event L occurs at time = VT— v 2 ® according 
to the out-bound observer and by eq. (177) at time tf = ^ ^j==g ) ^ according to the in-bound 
observer. By the way in which event A and B are defined we know that tp = tp and tf — tf. Finally, 
knowing that events A and B both occur at x E = 0, we see from eq. (176) that t E = (l — v 2 ) ^, and 
wc sec from eq. (177) that tf = (1 + v 2 ) -. Events occuring on earth between t E and tf are never 
simultaneous with either the out-going or in-coming observer on the rocket; there is consequently a 



sequence of events on earth during a time —t^ = 2vD; this precisely accounts for the inequality 
of the "twin paradox" noted above. 
From eqs. (166-170) it follows that 

t'* _ r '2 = f _ r 2 (17g) 

Eq. (178) is quite fundamental. We first note that it implies that the length of the four 
component vector x = (it, r) (where r is itself the spatial vector (x, y, z) and it is an imaginary 
fourth component) is unaltered by the boost transformations of eqs. (166-170), where the length of 
x is given by the usual dot product x.x. The boost transformation itself can then be taken to be a 
rotation in the it — x plane, albeit through an imaginary angle i6 where 

cos(i6) = 7, sin(i6) = ivy (179) 

It is more convenient though to define a real four component vector x^ = (t, r) and to combine the 
boost transformation and spatial rotations into the so-called "Lorentz" transformation. (If we also 
include translations in both time and space, we have the "Poincare" transformations.) The Lorentz 
transformations leave 

x 2 = i]^x v = t 2 - r 2 (180) 

unaltered. This four dimensional space with the magnitude of a vector being given by eq. (180) is 
called "Minkowski space" . The tensor (two component matrix) rj^ is diagonal with rj 00 = —rju = 
—V22 = ~ V33 = 1; it is the Minkowski space metric. (Indices associated with time are generally 
given the label "0".) It is convenient to define rf v to be the inverse of r/^ so that 

A general Lorentz transformation is given by 

a> = A^X 

Eq. (180) places a restriction on the transformation matrix A M v , 

If " is defined to be the inverse of A M v then it follows that 

A%A/ = £. (184) 
Eq. (182) now tells us that if <9 M = then 

d, = A\d' u ; (185) 



(181) 



(182) 



(183) 



eq. (184) converts eq. (185) to 

3, = A/ft,. (186) 

It is evident that d 2 = rj^d^dy is invariant under a Lorentz transformation. 

Vectors that transform like x^ in eq. (182) are called "contravariant" ; those that transform like 
<9 M in eq. (186) are called "covariant". (These types of vectors can be generalized to deal with 
transformations more general than the linear Lorentz transformations considered here.) A general 
tensor T^--f m which has m contravariant and n covariant indices transforms as follows 

From eqs. (183) and (184), it follows that 

A/ = iw^A", (188) 

and hence is a scalar. Furthermore, if T M is a contravariant vector, then T M = r]^ u T u is a covariant 
vector. 

It was noted by Thomas that two successive boosts that are not collinear do not in general 
constitute a boost, but rather are composed of a boost and a rotation. To see this, consider a 
boost in the z direction characterized by velocity (0,0, v) followed by a boost in the y — z plane 
characterized by velocity (0, 8v y , 5v z ) with 5vi being small. Two successive applications of eqs. (166) 
and (167) show that the resulting transformation is composed of a boost characterized by velocity 
(0, 5v y y/l — v 2 , v ) and a rotation about the x axis through a small angle 5uj ~ — (l — y/1 — v 2 ). 

We note that the Lorentz transformation of eqs. (166-170) do not immediately provide what 
an observer "sees" when observing an object, as it takes a finite amount of time for a light signal 
emanated by an object to reach the observer. If the signal were emitted at time r and received at 
time t by an observer at the origin, then since the velocity of light is constant to all observers 

c(t — t) — (distance travelled by the light signal 

in the observer's frame of reference). (189) 

If the object is moving with speed v along the x axis, being at the origin at time t — 0, then eq. 
(189) becomes 

c ( t _ ^L) = \ Xp \ , (190) 



v 

If x p > so that the object is receding from the observer, then by eq. (190) 



while if x p < so that the object is approaching the observer 



The apparent, or observed, velocity of the particle is v p = From eq. (191), we see that for a 
receding object, v p — > | as v — > c while from eq. (192) v p — > oo as v — > c for an approaching object. 

The notion of "constant accelerations" should also be considered in the context of special rel- 
ativity. It obviously cannot mean that the speed of an object increases indefinitely, as it does in 
Galilean relativity, as this would mean that the speed would eventually exceed c which contradicts 
special relativity. What we do mean by "constant acceleration" is that the acceleration of an object 
is the same according to all observers instantaneously at rest with respect to that object. We can 
now show that this criterion is met by an object following the trajectory 

x 2 - cH 2 = a 2 . (193) 

As will be discussed below, this trajectory results from application of a constant force. From eq. 
(193) it follows that 

*L = tl (194) 

dt (a 2 + c 2 t 2 ) 1 / 2 { } 

and 

d 2 x c 2 cH 2 

dt 2 (a 2 + c 2 t 2 y/ 2 (a 2 + c 2 t 2 yi 2 ' K ' 

From eq. (194), we see that the object is instantaneously at rest at time t = 0; by eq. (195) the 
acceleration of the object is c 2 /a at this time. But by eq. (178), the trajectory of this particle 
according to another observer moving with velocity v along the x axis is 

a* _ C H' 2 = a 2 . (196) 

This observer will hence see the object at rest at time t' = and at this time the acceleration will 
also be c 2 /a. Hence the "hyperbolic motion" described by eq. (193) is consistent with the criterion 
for constant acceleration in special relativity. We note that from eq. (194), || — > ±c and by eq. 
(195) -> as t -> ±oo. 

If such a uniformly accelerating particle is moving along the branch of the hyperbole 

4 - c 2 t 2 p = a 2 (197) 

which has x p > 0, then as r = ±^^x 2 — a 2 , eq. (189) becomes 

1 



/ - ( ±-^Jx 2 — a 2 



x p (198) 



where the positive sign in eq. (198) is for a receding particle and the negative sign is for an 
approaching particle. Solving eq. (198) we find that 



This is valid only for x p > and t > 0; if t < then no signal coming from the object can be 
detected by this observer who is located at the origin. The apparent velocity and acceleration of the 

2 2 

particle are | — and ^ respectively. This is consistent with eqs. (191) and (192), as the limit 
t — > oo corresponds to the object receding from the origin and the signal being emitted at a time 
t — > oo, while t — > + corresponds to the object approaching the origin and the signal being emitted 
at a time r — > — oo. At times r — > ±oo, the magnitude of the velocity of the object approaches c. 



1.7 The Point Particle 

We now consider the canonical formalism for a free particle. If a particle moves along a space- 
time trajectory specified by a vector rr^(A) with specified end points x M (Aj) and x M (A/), then the 
action for determining this trajectory can only be constructed out of the geometrical quantities that 
characterize this path. The simplest such quantity is just the path length of the trajectory. Since 
this length is given by ds 2 where 

ds 2 = rj^dx^dx" (200) 



our initial action is taken to be 



S = -m Vds 2 = -m / Vx 2 d\ (201) 



where m is a scale parameter (called the "mass"). Often yds 2 = dtV 1 — v 2 is called the "proper 
time" as it is an invariant that equals the time elapsed between two events that occur in the rest 
frame of the observer (ie, where v — 0). 

It is apparent from eq. (201) that the choice of A is arbitrary; one could replace A by A' = /(A) 
and leave the form of S in eq. (201) unaltered. This form of "gauge invariance" is associated with 
a first class constraint arising in the canonical treatment of S. To see this, we note that 

dyx 2 —mx a 

From eq. (202) we see that not all components of are independent as 

p 2 = m 2 ; (203) 

this is a first class primary constraint. There are no further constraints associated with this action. 
The canonical Hamiltonian of eq. (13) can easily be shown to vanish and thus by eq. (81) the 
Hamiltonian is entirely given by the constraint of eq. (203) 

H = u(p 2 -m 2 ). (204) 

The "time" derivative of a^(A) (that is, the derivative with respect to A) can be found from eq. 
(94) and thus we get 

x"(X) = 2up fM (205) 



as well as 

P, = 0. (206) 

From eq. (205) it is apparent that since u is arbitrary, x^ is arbitrary; this is a reflexion of the fact 
that A itself is arbitrary. 

The canonical formalism with the gauge choice x° = A can be pursued, but it is much easier to 
simply insert this choice of A into the action of eq. (201). This leads to 



S = —m Vl — v 2 dt = —m J 7 1 



dt (207) 



where we have set x^ = (t, r) and v = r. With this, we find that 



p = d tip^l = mv7 (208) 
ov 

and the Hamiltonian becomes 

H = pv - L = 7717. (209) 

(Inserting factors of c in eq. (209) when v = leads to the famous equation E = mc 2 .) From eqs. 
(207, 208) we see that 

H = Vp 2 + m 2 . (210) 

If we supplement the action of eq. (207) with a potential associated with a force F in the spatial 
direction x, we have the Lagrangian L = —my/1 — x 2 + Fx; eq. (3) then yields the equation of 
motion 

d 



= Fm (211) 

We see immediately that m'x = F when x — 0, that is, in the rest frame of the particle, we have 
Newton's equation F = ma. If x = x = when t — 0, then eq. (211) can be integrated to yield the 
trajectory given in eq. (193), 



t 2 -- 2 = -(j) 2 - (212) 



Thus we see that a constant force results in "hyperbolic motion" in Minkowski space. 

The action of eq. (201) disappears in the limit m — 0. Furthermore, the expressions for 
momentum and energy given by eqs. (208) and (209) respectively would vanish if v 2 were not to 
approach 1. In order to handle this massless limit, we introduce an auxiliary field e and consider 
the Lagrangian 

L = _I (^! + m 2 e \ . (213) 



2 V e 

The equation of motion for e is just 



which, when substituted back into eq. (213) just leads to the Lagrangian appearing in eq. (201). 
The advantage of introducing e is that it is possible to let m vanish in eq. (213) and retain realistic 
dynamics; it is apparent that in this massless limit the equation of motion for e implies that x 2 = 
so that massless particles move at the speed of light. If one were to make a canonical analysis of 
eq. (213), one encounters the primary constraint p e = for the momentum conjugate to e while for 
the momentum conjugate to x^ 1 is p^ = — — . The Hamiltonian becomes 

H = ^(-p 2 + m 2 ) (215) 

so that requiring that p e = generates the secondary first class constraint of eq. (203). The primary 
constraint p e = is also first class. With nl = p e and ir2 = \{p 2 — m 2 ), eqs. (118) and (119) lead 
to a gauge transformation generator 

G = X -e (-p 2 + m 2 ) + ip e . (216) 
This results in the infinitesmal changes 

Se = {G,e} = -e (217) 

and 

5x» = {G, x»} = ep» = . (218) 

e 

The transformations of eqs. (217,218) leave the action associated with the Lagrangian L of eq. 
(213) unaltered as L changes only by a total derivative, 



2d\ 



V m I e 

e 



(219) 



Arc length, used in eq. (201), is not the only quantity that can be employed to form an action for 
a particle moving along a trajectory x^{X)\ one can also use the extrinsic curvature of the trajectory. 
The action of eq. (201) is then modified to become 



= -/>( 



n, (^r) { (220) 



where again ds denotes the arc length ^dx^dx^. This action does not possess the reparameterization 
invariance of eq. (201). If we now express x M as a function of t = x° as in eq. (207), then 



d 2 x^ 

~d^^ V(l-v 2 ) 2 'l-v 2 ' (1-v) 2 



v.a a va.v 

+ T, T2 (221) 



where a = ^ and v = ^. The action itself is then 



S = - I' (m 7 - x + n (a 2 7 3 + (a. v) 2 7 5 ) ) dt. (222) 



Since this action depends on a = x, the equation of motion is given by eq. (11) and we have two 
canonical momenta 

dL 

n = — = -2^ ( 7 3 a + 7 5 va.v) (223) 
c/a 

and 

dL d dL . . 

p =^-*ar (224) 

From eq. (223) it follows that 

a = - (2/ry 3 ) _1 (tt - vtt.v) . (225) 

We also have 

v = x. (226) 

The Hamiltonian is then given by 

H (x, v; p, 7r) = p.x + 7r.v — L (x, x, x) , (227) 

which becomes 

H = p.v + m 7 " 1 - (2fi-f 3 y 1 (tt 2 - (v.tt) 2 ) . (228) 
The equations of motion for any quantity A are 

HA 

-^ = {AH} (229) 

where now the Poisson bracket is given by 

OA OB _ OA dB OA OB OA OB 
' <9x dp dp (9x dv dn dn dv 

With the Hamiltonian of eq. (228), it follows that L + S is conserved, where L, the "orbital" 
angular momentum, is given by Lj = eijkXjPk and S, the "spin" angular momentum, is given by 

Si — tijk'Vj'R'k ■ 

If circular motion is postulated, so that 

x = R (cos (tut) , sin (cut) , 0) (231) 

in a "rest" frame, then by eqs. (171-173) in a frame boosted with velocity u in a direction perpen- 
dicular to the plane of rotation, 

x = ^— vVl — u 2 sin(cut),vVl — u 2 cos(cut),uj (232) 

where v = cuR. Together, eqs. (228-230) lead to 

x = v (233) 

p = (234) 



v = - (/X7 3 ) 1 (tt - vtt.v) (235) 

and 

3 1 

7T = (7T 2 — (tT.v) ) V -7T.V7T — p + m7v. (236) 

2/ry ' /rp 

With eq. (232), eqs. (233-236) show that 

ii = (cos(wt), sin{ui), 0) (237) 
where a = 1 — v 2 and f3 = 1 — u 2 . Furthermore, it follows that 



provided that 



li 3 — a 

The numerical value of the Hamiltonian in eq. (228) is 



mu \/l — i , 

0,0, - -) 238 

WT^tf 1 + f 1 



9 2m a 2 B 

v = t, — ■ (239) 



m a/1 — v 2 



(240) 



L + S-10,0, '.' 2 ^7 1 (241) 



and the total angular momentum is 

- s = i n n 

V2 + t)Vl - v 2 

with L and S not being separately conserved. 

The above discussion shows that by incorporating the extrinsic curvature of the trajectory of 
a particle into its Lagrangian, some qualitative features of the intrinsic spin of a particle can be 
recovered. However, a more satisfactory classical realization of spin is obtained by making use of 
Grassmann coordinates. If xW an d ip^^X) are Grassmann variables, then the massless limit of eq. 
(213) can be generalized to 

I fx 2 ■ \ 

L = — - ( hiip.ip-ixx.ipj. (242) 

If a(A) is an arbitrary Grassmann function, then the Lagrangian of eq. (242) undergoes a change 

5L = 4-( d™lL) (243) 



d\ V 2e 

when 



5 X » = ia^, (244) 
5e = —iax, (245) 



6^ = -a ( — - ) (246) 



v e 2e 
and 

5x = 2a (247) 

and consequently the action itself is invariant. 

The momenta associated with e and x vanish; these are first class constraints which can be 
associated with the gauge conditions x — and e = — 1. The equations of motion then collapse 
down to the dynamical equations 

x = ip = (248) 

as well as the constraint equations 

x 2 = ±.i> = 0. (249) 

These follow from the Lagrangian 

L = - (x 2 — iip.ip) (250) 



2 

which shows that the momentum associated with ip^ is proportional to ip^ itself. This is a second 
class constraint whose Dirac bracket, determined from eqs. (91) and (151), imply that 

W,$ v y = iTf (251) 

An extension of eq. (242) to the case in which m 2 ^ requires an additional Grassmann field ip 5 . 
The Lagrangian becomes 



L=-i 
2 



— + m 2 e + i (^tp.ip — ip 5 ip 5 ^j — ix (^-^- + rmp b 



(252) 



Now the dynamical equations of motion are 

^ = ^ = ^ = 0. (253) 
The constraint formalism can be applied. The canonical momenta are 

— ( ^ - ^xVv) ( 254a ) 

(254b) 
(254c) 

(254d) 

(254e) 





dL 


-1 


Vp, = 


dx» 


e 




dL 


= 


Pe = 


de 




dL 




TV = 




2^ 




dL 


i 


7T5 = 






dip 5 


= ~2 




dL 


= 0, 


% = 


dx 



so that the canonical Hamiltonian is 

H c = - 



-(p 2 - m 2 ) + -x(p ■ i/> - mif> 5 ). 



(255) 



With the primary second class constraints of eqs. (254c, d) we obtain the Dirac Brackets (see eq. 
(91)) 



{A,B}* = {A,B} + i 



(256) 



from which it follows that 



{^5,^5}* = -i- 

The primary constraints of eqs. (254b,e) yield secondary constraints as 

{Pe,H c } = ^(p 2 -m 2 ) 
{tt x , H c } = --(p • i/> - mip 5 ) 



which are first class as 



{p • ip — mip 5 , p ■ ip — mifj 5 }* = i{p 2 — m 2 ); 



(257a) 
(257b) 



(258a) 
(258b) 

(259) 



there are no tertiary (third generation) constraints. 
The gauge generator of eq. (130) is now 



G = fi a (p a = B lPe + B(p 2 - m 2 ) + iF x n x + iF((p ■ - mi/> 5 ) 



(260) 



where (B,Bi) and (F, Fx) are Bosonic and Fermionic gauge functions respectively. Taking the total 
Hamiltonian to be 



H T = H C + U e ^, 



a 1 

= ~^(P 2 - ™ 2 ) + IjXkV ■ i> - m^ b ) + \ e p e + i\ x 7T x 

then eq. (131) becomes 

Bip e + B(p 2 — m 2 ) + iF 1 7c x + iF(p ■ ip — rrwp 5 ) 

+ {G, Ht}* — S\ a p e — i5\ x 7i x = 0. 



(261) 



(262) 



(We use the Dirac Bracket in eq. (262) which permits us to set the second class constraints equal 
to zero at the outset.) From eq. (262) we find that 

B 1 = 2B + iF X (263a) 

F 1 = -2iF. (263b) 

Taking G to be given by eqs. (260, 263), then the change 5 A = {A, G}* leads to 

<?x" = 2Bp» + iF^ = ^c" - ^xV^j + iFV (264a) 

be = 2B + iF X (264b) 

5x = 2F (264c) 

6^ = Fp» = ( x 1 " - -xV) (264d) 



2 V 2' 

5ip 5 = mF. (264e) 

These transformations have a group property determined by 

{G(B a ,F a ), G(B b , F b )}* = 2-^-(iF a F b )p e + (iF a F b )(p 2 — m 2 ). (265) 

We see that commuting two gauge transformations generated by G leads to a pure Bosonic gauge 
transformation with gauge generator B = iF a F b . 



